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Abstract 

We discuss a scenario in which the highest energy cosmic rays (CR's) and 

cosmological 7-ray bursts (GRB's) have a common origin. This scenario is 

consistent with the observed CR flux above 10 20 eV, provided that each burst 

produces similar energies in 7-rays and in CR's above 10 20 eV. Protons may 

be accelerated by Fermi's mechanism to energies ~ 10 20 eV in a dissipative, 

ultra-relativistic wind, with luminosity and Lorentz factor high enough to 

produce a GRB. For a homogeneous GRB distribution, this scenario predicts 

an isotropic, time-independent CR flux. 
PACS numbers: 98.70.Rz, 98.70.Sa 
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Three cosmic ray (CR) events with energies > 10 2O eV have recently been observed by 
the Yakutsk [1], the Fly's Eye [2] and the Agasa [3] experiments. These events pose two 
challenges to current models for the production of CR's. First, the high energies rule out 
most of the acceleration mechanisms so far discussed [4]. Second, since the distance traveled 
by such particles must be smaller than lOOMpc [5], their arrival directions are inconsistent 
with the position of any astrophysical object that is likely to produce high energy particles 
[6]. 

The origin of gamma-ray bursts (GRB's) is yet unknown. However, recent observations 
give increasing evidence that GRB's originate from cosmological sources [7] (for a differ- 
ent point of view see, e.g., [8]), and a variety of plausible cosmological sources have been 
suggested [9]. Whatever the ultimate source is, observations strongly suggest the following 
similar scenario for the production of GRB's [10]. The rapid rise time, ~ ms, observed in 
some bursts implies that the sources are compact, with a linear scale r ~ 10 7 cm. The 
high luminosity required for cosmological bursts, ~ 10 51 ergs _1 , then results in an initially 
optically thick (to pair creation) plasma of photons, electrons and positrons, which expands 
and accelerates to relativistic velocities [11]. This is true whether the energy is released 
instantaneously, i.e. over a time scale r /c, as a "fireball", or as a wind over a duration 
comparable with the entire burst duration (~seconds). In fact, the hardness of the observed 
photon spectra, which extends to ~ lOOMeV, implies that the 7-ray emitting region must 
be moving relativistically, with a Lorentz factor of order a few hundreds [12]. 

If the observed radiation is due to photons escaping the fireball/ wind as it becomes 
optically thin, two problems arise. First, the photon spectrum is quasi-thermal, in contrast 
with observations. Second, the plasma is expected to be "loaded" with baryons which may 
be injected with the radiation or present in the atmosphere surrounding the source. A small 
baryonic load, > 1O" 8 M , increases the optical depth (due to Thomson scattering) so that 
most of the radiation energy is converted to kinetic energy of the relativistically expanding 
baryons before the plasma becomes optically thin [13]. To overcome both problems Rees & 
Meszaros [14] proposed that the observed burst is produced once the kinetic energy of the 

2 



ultra-relativistic ejecta is re-randomized by some dissipation process at large radius, beyond 
the Thomson photosphere, and then radiated as 7-rays. Collision of the relativistic baryons 
with the inter-stellar medium [14], and internal collisions within the ejecta itself [15,16], were 
proposed as possible dissipation processes. 

In this letter we consider the possibility that high energy CR's and cosmological GRB's 
are related phenomena. We find that the observed flux of CR's beyond 10 2O eV is consistent 
with a scenario in which these particles are produced in GRB's provided that each burst 
produces similar energies in 7-rays and highest energy CR's. Adopting a dissipative ultra- 
relativistic wind model for the production of GRB's, we then demonstrate that protons may 
be accelerated in such a wind to energies of order 10 2O eV by Fermi's mechanism [17]. These 
two remarkable coincidences possibly suggest that GRB's and highest energy CR's have a 
common origin. 

GRB rate and CR flux. The high energy CR experiments observe a cone around the 
zenith of opening angle ~ 45° over a period of ~ lOyr. The rate of cosmological GRB events 
is v 1 ~ 3 X 10~ 8 Mpc~ 3 yr _1 [18]. The rate of events in a cone observed by CR experiments 
extending to a distance of ~ lOOMpc, which is the maximal distance likely to be traveled 
by a CR of 10 2O eV, is ~ 1 per 50yr. The rate in a cone extending to 50Mpc, the maximal 
distance traveled by the 3 X 10 2O eV CR observed by the Fly's Eye, is 8 times smaller. Thus, 
the probability of the experiments to observe a CR pulse from a GRB event over ~ lOyr 
period is small, unless the CR pulse is broadened in time, due to propagation through the 
inter-galactic medium (IGM), over a time scale > lOOyr. 

Time broadening is likely to occur due to the combined effects of deflection by random 
magnetic fields and energy dispersion of the particles. Consider a proton of energy E prop- 
agating through a magnetic field of strength B and correlation length A. As it travels a 
distance A, the proton is typically deflected by an angle a ~ X/Rl, where Rl = E / eB 
is the Larmor radius. The typical deflection angle for propagation over a distance d is 
6 S ~ (d/Xf/ 2 X/R L . This deflect ion results in a time delay, compared to propagation along 
a straight line, of order 
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r{E)^6ldlc^(^d)- c (1) 

((1) is valid for small 8 S ). The time delay is sensitive to the particle energy. Thus, a spread in 
particle energies AE ~ E results in a time broadening of the pulse over a time ~ t(E). Since 
the proton energy loss is a random process (collisions with the microwave background), there 
is a large dispersion in particle energies for a propagation distance over which significant 
energy loss occurs. For protons with energies > 10 2O eV the RMS energy spread is comparable 
to the average energy over propagation distances in the range 10 — lOOMpc [5]. Thus, the 
CR pulse would be spread over a time t(E). The field required to produce t > lOOyr is 

B > 2 x lO-^^o^Ai-o 172 G, (2) 

where E = 10 2 ° E 2 oeV , d = 100c?iooMpc and A = 10Ai O Mpc. The required field is not incon- 
sistent with our knowledge of IGM properties. Estimates of the average (over cosmological 
scales) IGM field are of order 3 X 10~ n [19], while disordered fields may be larger. A time 
broadening over t > lOOyr is therefore reasonable. 

The time variability and directional properties of CR flux produced by GRB events 
would depend on the distribution of GRB's in our vicinity and on the structure of the IGM 
magnetic fields. A detailed analysis of the flux properties are therefore beyond the scope 
of this letter. However, an order of magnitude estimate of the required CR production by 
a GRB may be obtained as follows. If GRB events are homogeneously distributed over 
scales of order ~ lOOMpc, and if t ^> lOOyr, then the CR flux due to GRB's would be 
approximately isotropic and time independent. The observed CR flux at energies beyond 
10 2O eV is ~ 3 X 10 17 m 2 s 1 sr 1 [20], corresponding (for isotropic time independent flux) 
to a number density ucr ~ 10~ 3O cm~ 3 of these particles. If this density is produced by 
GRB's then each GRB event should produce N = ucr/ 'v^tcr protons beyond 10 2O eV, 
where tcr ~ 3 X 10 8 yr is the life time of protons with energies > 10 2O eV [5]. The energy in 
protons beyond 10 2O eV, that should be released in each GRB in order to be consistent with 
the observed CR flux is therefore 
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E CR « 5 x 10 5O £ 2O erg. (3) 

Here, E = 10 2 ° E 2 oeV is the average energy of protons beyond 10 2O eV. Thus, Ecr is required 
to be similar to the energy emitted as 7-rays by a cosmological GRB, E 1 % 2 X 10 51 erg. 

We have so far assumed that the photon emission from GRB's is spherically symmetric. 
If the emission is concentrated into a solid angle 80,, then the rate of GRB events is higher 
by a factor 47r / 80, and the 7-ray energy produced by each event is smaller by the same 
factor. The above conclusions remain unchanged, except that Ecr in (3) is also reduced by 
a factor 47r/£fi. 

Fermi acceleration in dissipative wind models of GRB's. We consider a compact source 
producing a wind, characterized by an average luminosity and mass loss rate, denoted L and 
M = L/rjc 2 respectively. At small radius, the wind bulk Lorentz factor, 7, grows linearly 
with radius, until most of the wind energy is converted to kinetic energy and 7 saturates 
at 7 ~ r\. We assume that the wind is ultra-relativistic, 7 ^> 1, and that internal shocks 
in the ejecta reconvert a substantial part of the kinetic energy to internal energy at r ~ r^. 
Since 7 ^> 1, substantial dissipation of kinetic energy implies that the random motions in 
the wind rest frame are relativistic. The relativistic random motions are likely to give rise 
to a turbulent build up of magnetic fields, and therefore to second order Fermi acceleration 
of charged particles. We derive below the constraints that should be satisfied by the wind 
parameters in order to allow acceleration of protons to ~ 10 2O eV. 

The most restrictive requirement, which rules out the possibility of accelerating particles 
to energies ~ 10 2O eV in most astrophysical objects, is that the particle Larmor radius Rl 
should be smaller than the system size [4]. In our scenario we must apply a more stringent 
requirement, namely that Rl should be smaller than the largest scale R over which the 
magnetic field fluctuates, since otherwise the second order Fermi acceleration may not be 
efficient. We may estimate R as follows. The comoving time, i.e. the time measured in the 
wind rest frame, is t = r/'yc. Thus, at r = regions separated by a comoving distance 
larger than 7^/7 are causally disconnected, and the wind properties fluctuate over comoving 
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length scales up to R ~ ra/'y. We must therefore require Rl < ra/'y. A somewhat more 
stringent requirement is related to the wind expansion. Due to expansion the internal energy 
is decreasing and therefore available for CR acceleration (as well as for 7-ray production) 
only over a comoving time tj, ~ rj/'yc. The typical Fermi acceleration time is t a = fRi,jcf5 2 , 
where /3c is the Alfven velocity and / ~ few [21,4]. In our scenario (3 ~ 1 leading to 
the requirement /Rl < T^/'y. This condition sets a lower limit for the required comoving 
magnetic field strength. In obtaining this limit one should note, that due to the relativistic 
motion the energy of the accelerated particle is ~ 7 times smaller in the wind rest frame 
than in the observer frame. Thus, the Larmor radius is given by Rl = E/"feB, where E is 
the observed particle energy. Using this expression for Rl we obtain 

B>f— = 3x 10 4 /£2or^ 3 G, (4a) 
er d 

where r& = 10 13 rd,i3cm. 4a can be stated in a different form by comparing the required 
field to the equipartition field, i.e. to a field with comoving energy density similar to that 
associated with the random energy of the baryons. Substantial dissipation of kinetic energy 
implies that the contribution of the random motion energy density to the luminosity is large. 
In this case (4a) may be put in the form 

(^JyO-lSf^E^, (4b) 

where B ep . is the equipartition field, 7 = 3OO7300J and L = lO^i^iergs -1 . 

The accelerated proton energy is also limited by energy loss due to synchrotron radiation 
and interaction with the wind photons. The condition that the synchrotron loss time, 
t sy = (QirrripC 3 / (TTrnDE -1 B~ 2 , should be smaller than the acceleration time t a is 

B < 3 x lOV-Saoo^G. (5) 

This condition may be satisfied simultaneously with (4a) provided that the dissipation radius 
is large enough, i.e. 

r d > 10 12 / 2 73oo^2ocm. (6) 
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The high energy protons we are interested in lose energy in interaction with the wind 
photons mainly through collisions producing pions, for which ~ 10% of the energy is lost 
per collision [5]. The typical time for energy loss in this process is t 1 ~ 10/ra 7 o"c, where 
a % 10~ 29 cm 2 is the collision cross section for high energy protons and ra 7 is the photon 
number density. ra 7 is approximately given by L 1 = 47rr 2 c7ra 7 £ 7 , where e 7 is the observed 
photon energy, e 7 ~ MeV, and L 1 is the gamma-ray luminosity. The condition t 1 > t a gives 

B > 20/ J L 7 , 5lJ B 2 or d - 2 3 73o 2 oG. (7) 

Let us discuss the implications of (4-7) for a wind which should produce a GRB. The bulk 
Lorentz factor of the wind at must satisfy 7300 > 1. If the process converting the random 
energy to gamma-rays is not very inefficient, then L is a few times higher than the observed 
gamma-ray luminosity, L 1 % 10 51 ergs _1 . (4b) then implies that the magnetic field energy 
should be close to equipartition. This requirement is independent of the value of and 
determined only by the wind luminosity and Lorentz factor. The condition (5) is satisfied 
in this case provided that dissipation occurs at large enough radius, as required by (6) (this 
also ensures that (7) is satisfied). In the dissipative wind models for cosmological GRB's, 
proposed by Paczyhski & Xu [15] and by Rees & Meszaros [16], internal shocks occur due to 
time variability of the source. Variability that leads to fluctuations in the bulk Lorentz factor 
A7 ~ 7 over a time scale At, produces internal shocks which dissipate a substantial fraction 
of the wind kinetic energy at % "f 2 cAt. The flux emitted by the source may fluctuate 
on time scales ranging from the source dynamical time, td yn ~ r /c ~ 10~ 3 s, to the wind 
duration T ~ s, leading to dissipation radii in the range 3 X 10 12 7| oo cm to 3 X 10 15 7| oo cm. 
The dissipation radii in these models may therefore be large enough to satisfy (6) and allow 
the production of high energy protons. 

The dissipation radius of a GRB producing wind is limited by the observed duration, 
T ~ s, of the 7-ray pulse, < "f 2 cT. This condition may be satisfied simultaneously with 
(6) provided that the Lorentz factor satisfies 

7 > 40/ 1 / 2 £ 2 3 o /4 Tr 1/4 , (8) 
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where T = T\S. 

The estimates given above were obtained assuming that the wind is spherically symmet- 
ric. However, since the wind is ultra-relativistic, the behavior of a jet-like wind is similar 
to that of a spherical one, provided that the jet opening angle is larger than I/7. In this 
case, L should be understood as the luminosity that would have been produced by the wind 
if it were spherically symmetric. A possible complication that may arise for jet-like wind is 
that particles may escape the acceleration region transversally on a time scale shorter than 
t a (the condition Rl < ra/'y ensures that the radial escape time is larger than t a ). However, 
for jets with opening angles > I/7 the transversal size of the jet is larger than 7^/7, so that 
the transversal escape time is also larger than t a . 

In conclusion, we have found that the observed flux of CR's with energies > 10 2O eV 
is consistent with a scenario in which these particles are produced in cosmological GRB's 
provided that each burst produces similar energies in gamma-rays and CR's beyond 10 2O eV. 
We have also demonstrated that a dissipative ultra-relativistic wind, with luminosity and 
bulk Lorentz factor implied by GRB observations, satisfies the constraints necessary to allow 
the acceleration of protons to energies ~ 10 2O eV by second order Fermi acceleration. We 
emphasize that this energy estimate is independent of the wind model details, and depends 
only on the wind luminosity and average Lorentz factor, which are directly estimated from 
observations (assuming a cosmological origin for GRB's). 

Acceleration to ~ 10 2O eV is possible provided that the dissipation of kinetic energy 
occurs at radii > 10 13 cm [(6)], and that turbulent build up of magnetic fields during 
dissipation gives rise to fields which are close to equipartition [(4b)]. Although the details 
of the mechanism of turbulent build up of an equipartition field are not fully understood, it 
seems to operate in a variety of astrophysical systems. We note that a magnetic field of this 
strength may exist in the plasma prior to the onset of internal collisions if a substantial part 
of the wind luminosity is initially, i.e. at r ~ r , provided by magnetic field energy density, 
and if the field becomes predominantly transverse. The pre-existing field may suppress 
small scale turbulent motions [16]. In this case shocks coherent over a scale R ~ r^/7 would 
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exist and protons would be accelerated by first rather than second order Fermi mechanism. 
The constraints (4-7) are valid in this case too, therefore leaving the above conclusions 
unchanged. 

Future investigations of a scenario, in which the highest energy CR's are produced by 
cosmological GRB's, should aim at obtaining characteristic signatures of such a scenario, 
which may be compared with observations. A homogeneous distribution of GRB events over 
scales ~ lOOMpc would naturally lead to isotropic time independent CR flux. This predic- 
tion may be compared to observations once larger high energy CR detectors are operative. 
The CR energy spectrum, that should be produced by GRB's, is already constrained by ob- 
servations. The life time of protons with energies > 10 19 eV is larger than the Hubble time, 
~ 10 lo yr, which is ~ 30 times the life time of protons above 10 2O eV (~ 3 X 10 8 yr). Thus, if 
the number of protons produced above 10 19 eV in a GRB event is similar to the number of 
protons produced above 10 2O eV, then the ratio of the flux above 10 19 eV to the flux above 
10 2O eV would be ~ 30. The observed ratio is ~ 300. Thus, if the number of CR's produced 
per energy interval is a power law dN oc E~ n dE, as typically expected from Fermi accelera- 
tion [4], we must have n < 2. A more detailed analysis of the acceleration process, resulting 
in a prediction for the produced energy spectrum, may therefore be directly compared with 
observational constraints. Note, that if the power law spectrum extends from 10 21 eV down 
to 10 12 eV (the lower bound is due to the fact that the source bulk Lorentz factor is ~ 10 3 ), 
then n < 2 implies that the total energy of relativistic protons is < 5 X 10 51 erg [cf. (3)]. Since 
the observed gamma-ray energy is % 2 X 10 51 erg, we do not expect the energy associated 
with CR's to change significantly the constraint on the total GRB source energy. 

I thank J. W. Cronin for introducing me to the subject of highest energy cosmic rays, 
J. N. Bahcall and R. M. Kulsrud for invaluable comments, and F. Rasio and J. Goodman 
for helpful discussions. This research was partially supported by a W. M. Keck Foundation 
grant and NSF grant PHY 92-45317. 
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